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ABSTRACT 

Aims: We assess the sensitivity of void shapes to the nature of dark energy that was pointed out 
in recent studies. We investigate whether or not void shapes are useable as an observational 
probe in galaxy redshift surveys. We focus on the evolution of the mean void ellipticity and 
its underlying physical cause. Methods: We analyse the morphological properties of voids in 
five sets of cosmological N-body simulations, each with a different nature of dark energy. 
Comparing voids in the dark matter distribution to those in the halo population, we address 
the question of whether galaxy redshift surveys yield sufficiently accurate void morphologies. 
Voids are identified using the parameter free Watershed Void Finder. The effect of redshift 
distortions is investigated as well. Results: We confirm the statistically significant sensitivity 
of voids in the dark matter distribution. We identify the level of clustering as measured by 
as(z) as the main cause of differences in mean void shape (e). We find that in the halo and/or 
galaxy distribution it is practically unfeasible to distinguish at a statistically significant level 
between the various cosmologies due to the sparsity and spatial bias of the sample. 

Key words: Cosmology: theory - large-scale structure of universe - dark energy - methods: 
data analysis - numerical simulations 



1 INTRODUCTION 

One of today's greatest puzzles is that of the nature of the dark com- 
ponents of our universe: dark matter (DM) and dark energy (DE). 
According to concordance cosmology, these make up respectively 
21.7% and 73.8% of the total cosmic energy content ( [Komatsu] 
|et al.|20lT] ), a total of 95.5%. Yet, there are no compelling obser- 
vational indications for the precise nature of dark matter and dark 
energy ( [Komatsu et al.|2011[|Amanullah et al.|2010| ). 

While there is considerable agreement on some of the physi- 
cal properties of dark matter, the nature of dark energy remains a 
complete mystery. It is not even sure whether it really concerns an 
energy component to be associated to a new species in the universe 
or rather a modification of gravity itself, as in extensions of general 
relativity like scalar- tensor theories. 

It is far from trivial to constrain the nature of dark energy, due 
to the relatively weak imprint of the equation of state of dark en- 
ergy, in combination with sizeable observational errors (Fri eman,| 
|Turner &~H uterer 2008 ). Nonetheless, since the Supernova Cos- 
mology Project and the High-z Supernova Search Team first hinted 
to the accelerated expansion of the universe (Riess et al. 1998| |Perr| 
|mutter et al.||1999| ), a range of astrophysical probes have led to a 
considerable narrowing of the parameter space for the nature of 
dark energy. Weak gravitational lensing by the large scale matter 
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Figure 1. ACDM (left) and quintessence (right) cosmologies; dark matter 
density (grey) and galaxies (coloured dots). Differences are subtle; measur- 
ing them is the challenge of this paper. 
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distribution, baryonic acoustic oscillations (BAOs) and the redshift 
dependent number density of clusters of galaxies are well-known 
examples of dark energy probes. Based on these, a range of obser- 
vational and experimental programs have been initiated in order to 
constrain the nature of dark matter and dark energy. Notable exam- 
ples are BAO experiments like WiggleZ (Drinkwater et al. 2010), 
BOSS ( [Ross et al.|2010| ), the Dark Energy Survey and large weak 
lensing surveys like KiDS dde Jong et al.|2 012). Ambitious projects 
like ESAs Euclid Mission ( Laureijs ^et al.|201 1| ) and the Large Syn- 
optic Survey Telescope (LSST) will considerably sharpen and ex- 
tend our knowledge on the dark energy content of the universe. 
Most of these probes directly relate to the cosmic distance mea- 
sure. 

Probes like the integrated Sachs Wolfe effect (ISWWDupe 
et al.|20lT| and th e structure formation growth rate ( |Linder|2 006 ; 



Guzzo et al. 2008) offer alternative paths towards constraining dark 



energy. However, these probes are still ridden by considerable un- 
certainties. In preparation for the upcoming large dark energy ex- 
periments, it is therefore of great importance to find a wider range 
of reliable and independent probes for the determination of the 
dark energy equation of state. The combination of different and 
independent observations and measures of dark energy is crucial 
for breaking degeneracies and constraining the allowed parameter 
space. The recent study by Aman ullah et al.| ( |2010| > illustrates the 
potential of such an approach to further pin down dark energy and 
other cosmological parameters. 



1.1 Voids & cosmology 

Here, we are particularly interested in the imprint of dynamical 
dark energy ( |Wetterich|1988| >. 

Several recent studies have pointed out that cosmic voids not 
only represent a key constituent of the cosmic mass distribution, 
but that they are also one of the cleanest probes and measures of 
global cosmological parameters. Their structure, morphology and 
dynamics reflects the nature of dark energy (|Park&Lee|2007l[Lee1 



& Park 2009, Lavaux & Wandelt 2010, Biswas, Alizadeh & Wan- 



delt 2010 , Lavaux & Wandelt 201 1 , Shoji & Lee 2012), dark matter 
( [Hell wing, Jusz kiewicz & van de W eygaert 2010 , Li 201 1) and that 
of the possibly non-Gaussian nature of the primordial perturbation 
field ( Kamionkowski, Verde & Jimenez 2009). Lavaux & Wandelt 
( |2010| ) demonstrated the extreme sensitivity of the evolving ellip- 
ticity of voids to the equation of state of dark energy. [Biswas, AT] 
izadeh & Wandelt ( 2010) even quoted the possibility of improving 
the Dark Energy Task Force figure of merit by a factor of a hundred 
for future experiments like Euclid. 

Along a related and perhaps even more promising route, 
Lavaux & Wandelt ( 2011) demonstrated the potential of using an 
Alcock-Paczysnki test on the average shape of stacked voids (also 
see |Shoji & Lee|20"T2| . By stacking a large number of voids, one 
expects a spherical shape, so that size differences in radial and 
transverse direction can be directly related to the product of an- 
gular diameter distance and Hubble parameter. The claim is that 
for Euclid-like surveys stacked voids would outperform BAOs by 
an order of magnitude in accuracy ( Lavaux & Wandelt 2 011| ). 

Other studies have explored the possibility to use voids for in- 
ferring information on the amount and nature of dark matter. Void 
outflow velocities and the accompanying redshift distortions might 
be used to determine fi m , the mass density of the universe, and in- 
fer the amount of dark matter ( Martel & Wasserman 1990, Dekel 
|& Rees|1994]|Ryden & M elott 1996). However, void outflow ve- 
locities are difficult to measure, while void redshift distortions are 



restricted because of the naturally limited void outflow. Nonethe- 
less, recent advances in the study of cosmic flows, have shown the 
prominent and clearly recognizable dynamical role of voids in the 
nearby universe ( |Courtois et al.|2012[ |Tully et al.(2008) . Related 
observations have been made with respect to the sensitivity of voids 
to the identification of dark matter. The emptiness of voids would 
be a direct measure of the strength and screening length of the class 
of long-range scalar-interacting dark matter models that have been 
forwarded as a possible means of explaining several deficiencies 
of the concordance ACDM model ( Farrar & Peebles 2004, Gub- 
|ser & Peeblesl|2004l |Nusser, Gubser & Peebles 2005 , Peebl es~& 
Nusser 2010, Hellwing, Juszkiewicz & van de Weygaert 20101 ID" 



201 1| >. Most prominent amongst these deficiencies is the observed 
dearth of dwarf galaxies in nearby voids ( Peebles 2001). Given the 
extreme environment of voids, probing the tail of the cosmic den- 
sity and halo mass distribution (Platen, van de Weygaert & Jones 
|2012|), they form a natural resort for exploring the imprint of pos- 
sible modifications of general relativity, such as f(R) gravity < [0| 
[20TT| > and MOND/TeVeS dLlinares|20lf] ». 



1.2 Voids 

Voids form a prominent aspect of the megaparsec distribution of 
galaxies a nd matter ([Chin carini & Rood|1975[|Gregory & Th omp- 
son! 19781 [Einasto, Joeveer & Saar|198o[|Kirshner et al.|19811|KiT 



shner et al. 1987^ Ide Lapparent, Geller & Huchrall 19861 Ivan de 
Weyg aert|1991||Colless et al.|2003||Tegmark et al.|20041 |Furlan 
etto & Piran 2006; Huchra et al. 2012). They are enormous regions 
with sizes in the range of 20 — 50 /i _1 Mpc that are practically 
devoid of any galaxy. They are usually roundish in shape and oc- 
cupy the major share of space in the universe (see van de Weygaert 
|& Platen|2011| for a recent review). Surrounded by elongated fila- 
ments, sheetlike walls and dense compact clusters, they weave the 
salient weblike pattern of galaxies and matter pervading the observ- 
able universe. 

Voids in the galaxy distribution account for about 95% of the 
total volume ( Joeveer, Einasto & Tago 1978 Kauffmann & Fairall 
|E1-Ad, Piran & da Costa|1996||El-Ad & Piran|1997| Roj as 



, 2005 , Pan et al. 2012). The typical sizes of voids in the galaxy 



distribution depend on the galaxy population used to define the 
voids. Voids defined by galaxies brighter than a typical L* galaxy 
tend to have diameters of order 10 — 20 /i _1 Mpc, but voids associ- 
ated with rare luminous galaxies can be considerably larger; diam- 
eters in the range of 20 — 50 h~ 1 Mpc are not uncommon |Hoyle| 
& Vogeley 2002, Plionis & Basilakos 2002). Even larger voids can 
be recognized in the distribution of clusters ( Bahcall 1988 , Einasto 
|etal.|1994||200l) . 



Evolving out of primordial underdensities, voids become in- 
creasingly isotropic objects (Icke 1984) with a "bucket- shape" den- 
sity profile whose density in the centre has a characteristic value of 
5 « —0.8 (Sheth & van de Weygaert 2004). To a first approxima- 
tion, (isolated) spherical underdensities will become more spherical 
as they expand ( Icke 1984 , van de Weygaert & van Kampen 1993). 
In reality, voids will never reach perfect sphericity. Their flattening 
is a result of large scale dynamical and environmental influences 
( Platen, van de Weygaert & Jones 2008). They will encounter sur- 
rounding structures such as overdense filaments or walls. More- 
over, they retain an uncommonly large sensitivity to the dynamical 
influence of their large scale environment. In most situations this 
remains the dominant factor, to the extent that voids are found to 
become more anisotropic as time proceeds. Under realistic circum- 
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stances, the evolution of voids appears to reverse the simple spher- 
ical trend expected for isolated voids (Icke |1984| . 

The sensitivity of voids to global cosmological parameters is 
a result of their unique dynamical status. On the one hand, the dy- 
namical evolution of voids is rather straightforward, in that they 
tend to evolve into expanding, extended, uniform, and underdense 
regions with a distinct bucket- shaped profile ( Shet h & van de Wey^| 
gaert 2004 ). On the other hand, they are distinctly nonlinear ob- 
jects that mark the transition scale between linear and nonlinear 
evolution (Sahni, Sathyaprakah & Shandarin 1994). As such, their 
structure and morphology reflects and magnifies cosmological dif- 
ferences present in the primordial universe. Also, unlike the ma- 
jority of evolving overdensities in the form of dark matter haloes 
and galaxies, their evolution retains the dominant influence of the 
inhomogeneous cosmological surroundings (Plate n, van de Wey^] 
gaert & Jones 2008). All these factors conspire to make voids into 
important and optimal cosmological sources of information. 



While moderate density regions (5 ~ 0— a few) appear to have a 
rather low level of bias, a range of studies have indicated a strong 
bias of galaxies in and around the lowest density regions. Peebles 
(2002) strongly emphasized the fact that dwarf galaxies seem to 
avoid the void regions, contrary to the expectations and predictions 
of most galaxy formation theories. This void phenomenon is one 
manifestation of what appears to be a rather puzzling situation with 
respect to the galaxy distribution in voids: different semi- analytical 
schemes of galaxy formation predict either bias or anti-bias ( Platen 
|2009l [Platen, vg n de Weygaert & Jones|2012| ). 

In figure [T]the galaxies/haloes (points) are overlaid on the full 
density field, illustrating that the spatial distribution of the two 
don't match. This is an especially poignant example of the fact 
that void substructure is lost on small scales. The large void of 
~ 30 h~ 1 Mpc at the right of the figure is traced only on the outside 
by galaxy-like "haloes" (more on this in §3.2| ). 



1.3 Observing voids 

The use of voids as a cosmological probe involves at least three 
major complicating factors. 

A first point of attention is the very definition of a void. This 
involves the need for a clear and unambiguous detection and de- 
lineation of the near empty void regions. A range of studies have 
forwarded a large and diverse number of techniques and methods 
to accomplish this. For a complete overview of these algorithms 
we refer to the comparison study by Colber g et al.| < |2008| > and to 
the study of Lavaux & Wandelt (2010). We will use the Watershed 
Void Finder, developed by our group ( [Platen, van de W eygaert & 
Jones 2007). This method allows a parameter-free determination of 
the size and shape of voids in the matter and galaxy distribution. 
We will discuss some of the relevant technical details in £|4] 

A second issue concerns the fact that the relevant observations 
of cosmic structure consist almost exclusively of galaxy redshift 
surveys. The use of redshifts for distance estimation distorts the 
observed shapes of structures (Shoji & Lee 2012). This is a conse- 
quence of the peculiar velocities in and around these objects. In the 
case of voids, matter and galaxies are flowing out of the density de- 
pressions as a result of the lower than average gravitational attrac- 
tion. As a result, with respect to the Hubble flow we see galaxies 
at the frontside moving towards us and those at the backside mov- 
ing away from us, stretching the void along the radial direction. A 
number of studies in the nineties proposed this as a possible means 
of extracting global cosmological information from voids (e.g|Ry-| 



den & Melott 1996). Recent studies confirm this idea ( Percival & 



White 2009 



[Jennin gs, Baugh & Pascoli 2011| ). 



A third and major complicating factor for the study of voids 
in observations is the fact that we have to infer their properties 
from the diluted and possibly strongly biased population of galax- 
ies. Galaxy redshift surveys are necessarily limited in their spa- 
tial resolution of cosmic structure. Even in the Sloan Digital Sky 
Survey (SDSS), structure is delineated by relatively bright galax- 
ies. The spatial void distribution defines an intrinsically hierarchi- 
cal complex of voids, in which a void consists of an assembly of 



subvoids ( Sahni, Sathyaprakah & Shandarin 1994, Sheth & van de 
Weygaert|2004| |Aragon-Calvo et al.|2010[|Aragon-Calvo & Sza- 



lay 2012). Finer substructure can only be observed when sampling 



ever fainter objects. Hence, the brighter galaxies of SDSS only 
trace the rough outline of larger voids, with sizes in the order of 
10 — 15 /i _1 Mpc. Even more complicating is the issue of the pos- 
sible bias of galaxies with respect to the dark matter distribution. 



1.4 Voids in the halo distribution 

In this study, we address the issue of whether it is feasible to infer 
the nature of dynamical dark energy from the evolving population 
of voids in the observational context of the voids being sampled by 
a dilute population of discrete objects, like dark matter haloes and 
galaxies. 

In this, we follow the strategy of identifying voids only and 
solely on the basis of the spatial distribution of the sampled ob- 
jects. From the distribution of galaxies or haloes, we outline voids 
with the help of the Watershed Void Finder (WVF). While under 
ideal circumstances, one might consider to use knowledge of the 
dynamical evolution of voids, we chose not to do so given the un- 
certainties about the bias of void galaxies: it would amplify any 
difference between the sampled halo/galaxy distribution and the 
underlying density field. Our plan is to follow the evolution of the 
shape of voids inferred from the spatial dark matter and dark halo 
distribution at each epoch and to compare the results obtained for 
different dynamical dark energy models. The void population in the 
dark matter distribution of the different dark energy cosmologies is 
expected to reflect the sensitivity to the nature of dark energy. How 
accurately the void population in the dark halo distribution man- 
ages to follow this sensitivity, is one of the main questions to be 
answered by this study. 

Also, we investigate the origin of the dependence of voids to 
the nature of dark energy. We presume that the structure of voids at 
any cosmic epoch is a reflection of the stage of structure develop- 
ment in the corresponding dark energy cosmology. This suggests 
a strong relation between the amplitude of the density and veloc- 
ity perturbations in the mass distribution and the structure, shape 
and profile of the voids. We therefore assess in how far the void 
shape evolution in different dark energy scenarios can be ascribed 
to as, the root mean square density fluctuation within a sphere of 
radius 8 /i _1 Mpc, and whether there are other possible factors of 
influence. 



1.5 Outline 

The outline of this paper is as follows. In ^2] we will treat cosmo- 
logical models, four of which contain a time dependent dark energy 
component and discuss cosmic voids and their theoretical use in 
discerning between cosmological models. We will give a descrip- 
tion of the data in ^3] Void identification is invaluable to our anal- 
ysis and we will elaborate on this in ^4] In ^5] we describe the void 
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Figure 2. Equation of state parameter w versus redshift. Shows the evolu- 
tion of the dark energy equation of state parameter w = P/p. The different 
lines indicate different dark energy models. 



populations in the simulations of the various dark energy scenar- 
ios, along with their size and shape characteristics. Subsequently, 
in ^6] we examine the evolution of the void shapes, and try to as- 
sess under which conditions these can be used as probes, given the 
discrete, diluted and biased halo and galaxy samples on which the 
analysis is based. We finish our analysis in ^7] where we touch on 
the issue of clustering dependence. We conclude and discuss our 
findings in ^8] 



2 DARK ENERGY & VOIDS 
2.1 Dark energy cosmology 

The precise nature of dark energy will be a decisive factor for the 
fate of the universe. Our reference point in comparing different cos- 
mological models is that of "concordance" cosmology or ACDM 
cosmology. The geometry of the universe is flat and its matter con- 
tent is dominated by a species of cold dark matter, while baryonic 
matter represents a smaller proportion. Its dynamics is currently 
dominated by the cosmological constant A. 

In this paper, we are going to investigate a range of possible 
dynamical dark energy alternatives to the cosmological constant. 
The nature of dark energy is often specified in terms of the equation 
of state, w = P/p. The first order parametrization of its equation 
of state, 



w(a) = w + w a (l - a) . 



(1) 



is often used to compare different dark energy models (Chevallier 
|& Polarski|200T| |Linder|2003"] >. The impact of the nature of the - 
possibly dynamic and evolving - dark energy on the evolution of 
the universe, follows directly from the Friedmann equation, 



H_ 



exp 



if! 

J a Q 



+ Wj{a!) 



da 



(2) 



In this expression the Hubble parameter is H = a/a, a(t) is the 
(normalised) cosmological scale factor and the dot means deriva- 
tive with respect to cosmic time. The relative energy densities of the 
cosmological fluids - such as baryons, dark matter, radiation and 
dark energy - at the current cosmological epoch is given by floj. 
For a flat universe, we evidently have J2j — ^o,j — 1- 



2.2 Dynamical dark energy 

There are two enigmas if dark energy is in fact to be ascribed to 
a cosmological constant A. The first riddle concerns the realisation 
that it is relatively unlikely for the measured densities of matter and 
dark energy to be almost equal at the present time, because - ac- 
cording to a theory with a cosmological constant - they must have 
been very different in the past and they will be very different in 
the future. This is aptly called the coincidence problem, and con- 
cerns the question why the density of matter and dark energy are 
of comparable magnitude at the present epoch. The second puz- 
zle concerns the embarrassingly small value of the cosmological 
constant. Particle physics tells us that there should be a vacuum en- 
ergy, and that this would have the effect of a cosmological constant. 
However, the predicted value of the energy density of the vacuum 
is about 120 orders of magnitude higher than what is measured for 
Qa. If the interpretation of the cosmological constant as a form of 
vacuum energy would be right, this would call for a fine tuning of 
the initial conditions over an unconceivable large dynamic range. 
This is known as the fine-tuning problem. 

Dynamical dark energy has been proposed in order to find a 
possible and more natural solution ( Wetterich 1988, Ratra & Pee- 
|bles| [l988). The most interesting cases are obtained in the pres- 
ence of attractor solutions (Liddle & Scherrer 1999), in which the 
dark energy scalar field follows a trajectory which is the same for 
a wide range of initial conditions of the scalar field and its deriva- 
tive. Depending on the potential, dark energy can track the back- 
ground component and then dominate over it at relative recent times 
( |Copeland, Liddle & Wands|1998| ). 

In this work we investigate four cosmologies with different 
time dependent dark energy models. The first two models are basic 
quintessence models: the Ratra-Peebles (RP) and SUGRA models. 
These models involve a scalar field that does not interact with the 
other cosmic species, except minimally through the overall cosmic 
expansion. The difference between RP and SUGRA is in the po- 
tential V(4>) that describes the dynamics of the scalar field (j). The 
other two models are extended quintessence (EQ) models. Their 
Lagrangian contains a term which is responsible for a non-minimal 
coupling of the scalar field to gravity. Effectively, this introduces 
a "fifth force" between matter particles which modifies their inter- 
action ([Baccigalupi, Matarrese & Perrotta||2000 



Perrotta & Bac- 



cigalupi 2002, Pettorino, Baccigalupi & Mangano 2005 , Pettorino, 
Baccigalupi & Perrotta 2005 , Pettorino & Baccigalupi 2008). Ef- 
fectively, the extended quintessence models are scalar- tensor the- 
ories of gravity ([Hw ang 1990a b, Faraoni 2000, Bois seau et al.| 



|2000HEsposito-Farese & Polarskii2001[|Riazuelo & Uzan|2002| ). 
The coupling constant can be either positive or negative, with dif- 
ferent effects. We model both possibilities; the model with negative 
constant is called EQn and the one with positive constant is called 
EQp. In addition to these four cosmologies, we also investigate the 
standard ACDM cosmology. 

The differences between the five dark energy models can be 
appreciated on the basis of the evolution of their equation of state 
parameter, w(z), which is shown in figure [5] The largest differ- 
ences, over the full redshift span, are those between ACDM and 
SUGRA. At z = 2, the difference is no less than 0.4. At z = 0.5 
it has decreased to 0.2. The technical details on the models are pre- 
sented in appendix [A] and in |De Boni et al.| ( [2011| ). 

Minimally coupled quintessence models, even including a 
tracking behaviour, still suffer from the problem of fine tuning 
of the initial conditions, because the observational bounds on the 
dark energy equation of state are increasingly converging towards 
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a value of wde very close to —1. The closer iude is to the cos- 
mological constant value, the smaller the range of allowed initial 
values for the corresponding scalar field has to be, since the dy- 
namics of such a field is extremely constrained by the flatness of 
the potential in which the field evolves (Matarrese, Baccigalupi & 
|Perrotta|2004l|Pettorino, Baccigalupi & Perrotta|2005| ). This is also 
the reason why people have been investigating further alternatives, 
often requiring extensions of general relativity which modify the 
gravitational attraction felt by matter particles. Two popular ways 
of proceeding are either modifying the coupling with gravity itself, 
in the Jordan frame, as in scalar-tensor, EQ, theories or to change 
the coupling of dark matter particles only, directly in the Einstein 
frame ( Amendola 2000, Pettorino & Baccigalupi 2008 and refer- 
ences therein). 

In this paper we follow the first path. In EQ cosmologies, the 
non-minimal coupling to gravity induces an "R-boost" mechanism 
responsible for early, enhanced scalar field dynamics, by virtue of 
which the residual imprint of a wide set of initial field values is 
cancelled out. These models, therefore, 'extend' the attractor so- 
lution behaviour of quintessence fields to scalar-tensor cosmolo- 
gies. However, we note here that even in these cases a flat potential 
is still required in order to get a reasonable equation of state to- 
day. Also, scalar-tensor cosmologies (EQ and f(R)) involve a cou- 
pling to baryons too, which therefore require the presence of some 
'chameleon' -like mechanism that protects the mass of the scalar 
field in high density regions. Nevertheless, this set of models has 
interesting solutions (attractors) and is simple to implement, there- 
fore being a good candidate to test differences with quintessence 
models through N-body simulations. 



2.3 Void shapes & dark energy 



Following the earlier suggestions by |Park & Lee| < [2007| > ; |Lee & 
Park (2009), recent studies by Wandelt and collaborators JCavaux 
& Wandelt|20T0l|Biswas, Alizadeh & Wandelt|2010l > showed that 
voids may be used as precision probes of dark energy. 

The sensitivity of voids to dark energy is a result of the way 
in which the dark energy equation of state affects the dynamical 
evolution of voids via its imprint on the large scale tidal force fields. 
Their influence remains important during the entire evolution of 
voids. The tidal forces evoke a significant anisotropic effect in the 
development of the voids, even sometimes causing their complete 
collapse ( Sheth & van de Weygaert 2004). 

As a result, the elliptical shape parameters that describe the 
flattening or elongation of a void are expected to be intimately con- 
nected to the local tidal tensor (see equations|4]&j5]andjBond, Kof- 
man & Pogosyan 1996 ; Park & Lee 2007 ; Platen, van de Weygaert 



& Jon es 2008). It relates their shapes directly to the surrounding in 
homogeneous cosmic matter distribution responsible for the gravi 
tational tidal field. In turn, the evolution of the tidal forces are di- 
rectly dependent on the nature of the dark energy content of the 
universe. This offers a path towards measuring the cosmological 
parameters. 

The dependence of the void's sphericity s and oblateness p 
(equation [15} on the tidal tensor Tij, defined as the traceless com- 
ponent of the second derivative of the gravitational potential 0, 



rp 1 2 



(3) 



values of the tidal tensor, Ai > A2 > A3: 

1 + (6 V -2)s 2 +p 2 



Ai(p,s) 



A 2 (p, s) 



p 2 + s 2 + 1 



1 + (5 V - 2)p 2 + s 2 
p 2 + s 2 + 1 



(4) 



(5) 



where S v = Y^h=i We can infer the dependence of the elliptic- 
ity distribution of voids on the cosmological model. The probability 
density distribution f(s) for the sphericity s = 1 — e, with e the 
ellipticity, is as follows: 



f(l-€]z) = f(8]Z,R L ) 

3375^2 



I 



J V\p,s\S = S v ' 1 a(z,R L )]dp 
5^ + 15^(Ai + A 2 



exp 



x exp 



Vld^a 5 (z,R L ) 
15(A?+AiA 2 + Ai^ 



2a 2 (z,R L ) 



2a 2 (z,R L ) 
(2Ai + A 2 - 5 V ) 



x (Ai - A 2 )(Ai + 2A 2 -S v ) 



A(5 V - 3) 2 ps 
(p 2 + v 2 + l) 3 



(6) 



This distribution is sensitive to changes in the cosmological 
parameters ( [Park & Lee||2007| ) through ct(z,Rl), the linear rms 
fluctuation of the matter density field smoothed on a scale of Rl at 
redshift z, defined as: 



(z,R L ) = D [ 



> / \ f k dk 

' {z) i 2* 



P(k)W 2 {kR L )d\nk, (7) 



where D(z) is the linear growth factor, W{kRi) is the top-hat 
window function and P(k) the linear power spectrum. The filter 
scale Rl is directly related to the void size in Lagrangian space, 

Equation [6]implies that the mean ellipticity of voids decreases 
with redshift z, where the mean ellipticity (e) is defined as 



/ 



e /(e, z)de . 



(8) 



can be directly inferred from their relation to the (ordered) eigen- 



More importantly, the model's rate of ellipticity decrease is sen- 
sitive to changes in the cosmological parameters. The redshift de- 
pendence of the mean ellipticity can, thus, be used to discriminate 
between different values of w a ( |Lee & Pa rk 2009 ; Lava ux & Wan-| 
|delt|2010) . 

These results support the impression that voids are a promis- 
ing probe of the nature of dark energy. Indeed, as Biswas, Alizadeh 
& Wandelt ( 2010 ) showed, we can improve upon our constraints on 
cosmological parameters by including void shape data, e.g. from 
Euclid. One of the central tests of this study is to see whether in- 
deed we can use this probe in our simulations. 



3 SIMULATIONS & SAMPLES 

Our analysis is based on numerical N-body simulations of DM par- 
ticles in five different cosmological backgrounds, four of which in- 
clude evolving dark energy. We use version 3 of the GADGET code 
(Springel 2005), which has the ability to specify the mode of evo- 
lution of dark energy. This is achieved through an extended dark 
energy implementation, described in Dolag et al. (2004). 

We simulate boxes with periodic boundary conditions. The 
physical linear size of a box is 300 h~ 1 Mpc. They contain DM 
particles only. The initial conditions at z = 60 are set up using 
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Table 1. Dark energy model parameters used in the simulations. For a de- 
scription of the models see ^A] a is the power-law slope of the quintessence 
potential, £ is the coupling constant for the extended quintessence models 
and wj bd is given by equation [A9] wq is fixed at z = 0. w a is determined 
using a x 2 fit (see equation[TJ and as is given at z = (<js is the same for 
all models at z ~ 1089). 
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Figure 3. as as a function of redshift. The standard deviation of the values 
of the density field on a scale of 8 h~ 1 Mpc, as, is a measure of the amount 
of clustering in the field. We show the values of as for the ACDM and four 
quintessence models as a function of redshift. This is an indicator of the 
evolution of structure in the different cosmologies. 



the Zel'dovich approximation (Zel'dovich 1970). The general cos- 
mological parameters adopted for the simulations are the WMAP 
3-year data values: Q m = 0.268, Q A = 0.732, Q b = 0.044, 
h = 0.704, a 8 = 0.776 and n = 0.947. We use two different mass 
resolutions (256 3 and 768 3 particles) for our simulations and build 
DM halo catalogs from the 768 3 sets. 

We simulate five different dark energy models ( §2.2| ). The spe- 
cific model parameters are summarised in table [T] These parame- 
ters are consistent with current observ ational constraints dAcqua-| 
|viva et al.|2005l|Amanullah et al.|2010||Komatsu et al.|2011) . The 
tables of w(a), needed for the extended dark energy implementa- 
tion, are calculated for the different models (see figure [2|. 

The dark energy models are normalised at the CMB using the 
relation 



0"8,DE — CT 8 



^acdm(^cmb) 
^de(^cmb) 



(9) 



where we assume zcmb = 1089 and D is the linear growth fac- 
tor, which is dependent on the dark energy model through H. This 
rescaling will cause differences in the amount of clustering, char- 
acterised by the normalisation parameter as. In figure [3] we show 
the resulting measured values of as for five different models (taken 
from our own low resolution data, see below). It is clear that the 
dark energy models have their influence on structure formation. In 
^5] and ^6] we will see that its impact is significant. We describe 
these samples in more detail below. 



3.1 DM particles 

For each one of the five different dark energy models ( §2.2[ ), we 
use a high resolution cosmological N-body simulation of 768 3 DM 
particles. Each of these simulations has exactly the same initial con- 
ditions. The particles have masses of 0.443 x 10 10 /i _1 M©. For the 
ACDM model we have snapshots at z = 0.1, z = 0.25, z = 0.51, 
z — 1.00 and z — 2.04, while for the other models we have a high 
resolution snapshot at z — 0. More information on these data can 
be found in |De Boni etaLl(2011) . 

For each of the five DE models we also generated and com- 
puted eight low-resolution simulations, each with 256 3 DM par- 
ticles. The physical parameters used for the initial conditions of 
these sets are exactly the same as those of the high resolution sim- 
ulations. The lower resolution simulations are used to investigate 
the time evolution of the dark energy models and to investigate the 
influence of mass resolution on the statistics used in our analyses 
(^6j). By using different random realisation for the initial condi- 
tions, these simulations are also essential for assessing and ruling 
out possible cosmic variance effects ( §4.3| ). 

In addition, we obtain random particle samples from the high 
resolution simulation sets, each sample containing the same num- 
ber of particles as the number of dark matter haloes. These samples 
are used to evaluate the effect of biasing of the halo population (^6|. 



3.2 DM haloes 

For every high resolution simulation box we use the SUBFIND al- 
gorithm of Springel et al. (2001) to find the gravitationally bound 
haloes (groups of adjacent particles, representing concentrated 
clumps of DM) that can colloquially be identified as galaxy haloes. 

The SUBFIND haloes (also called subhaloes, a term we will 
not use) trace the general structures present in the field. Figures [T] 
and|4]show that haloes are found in clusters, filaments and also the 
more tenuous walls. The voids remain largely empty. Moreover, 
careful observation shows that the more massive haloes are found 
in the higher density cluster nodes and most pronounced filaments. 
This reflects the fact that the mass function of haloes is dependent 
on the large scale environment: the mass function of void haloes is 
shifted to low mass objects (see e.g. |Aragon-Calvo et al.|2007) . 

Most important for our study is the observation that while the 
haloes do trace the outline of the most prominent large scale fea- 
tures, the detailed substructure is lost in the halo distribution. This 
is immediately obvious when comparing the spatial halo distribu- 
tion in figure[4]with the underling dark matter distribution. We may 
therefore expect to find a substantially different void size - and 
probably also void shape - distribution than in the intricate large 
scale structure visible in the dark matter distribution. 



3.2.1 Halo data 

We produce halo catalogues of all the dark energy simulations at 
redshifts of z = 0, z = 0.1, z 0.25, z 0.51, z 1.00, 
z = 2.04, z = 2.98 and z = 3.80. Because the halo finder finds 
significantly less haloes at higher redshifts, the corresponding halo 
samples are unsuitable for comparison with the lower redshift halo 
samples. Some basic data on the halo samples in the five different 
cosmological simulations are listed in table [2] 

The position of a halo is simply defined as its centre of gravity, 
which is the mean position of its N constituent particles, given each 
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Figure 4. Halo distribution. In the left panel we show the distribution of haloes (gold) and groups of haloes (red) in a slice of thickness 0.25 fo -1 Mpc of the 
ACDM simulation at z = 0. The density is shown in grayscale; the corresponding values are indicated by the color bar to the right. In the right zoom-in panel 
we additionally split the haloes by mass; red dots have M < 1O 12 M0, blue open diamonds have 1O 12 M0 < M < 1O 13 M0 and the two green filled circles 
are even more massive; the centre one has a mass of 2.4 • 1O 13 M0, the lower right one has 3.1 • 1O 14 M0. The orange circles indicate the location of haloes 
and groups. The radii of the circles are proportional to their mass. 
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Table 2. Characteristics of the halo catalogues obtained from the DE sim- 
ulations, at z = and z = 1. From left to right, we show the number of 
haloes A^h a i oes , the mean halo mass (m) , the standard deviation of the halo 
masses cr m and the largest mass in that snapshot m max ■ Masses are in units 
of 10 10 M o . 



has the same mass: 



1 N 



N 



(10) 



where x\ are the particle positions and N is the number of particles 
in the halo. The velocity is calculated analogously. Halo mass is 
simply the sum of the particle masses. 

In figure [5] we show some further characteristics of the halo 
catalogues as a function of redshift. The upper left panel shows 
the total number of haloes found at the given redshift. We observe 
that up until about z ~ 0.5 the number of haloes is largest in the 
ACDM cosmology. After this, the quintessence models take over. 
This hints at a lower rate of halo mergers, the main driving force 
behind hierarchical structure evolution. Indeed, it is mainly due to 
the lowest mass haloes, meaning that these have not been able to 
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Figure 5. Characteristics of the halo catalogues as a function of redshift. 
Different lines show catalogues obtained from the different dark energy 
simulations. Top left: number of haloes A^ioes • Top right: the mean halo 
mass (m). Bottom left: the standard deviation of the halo masses cr m . Bot- 
tom right: the largest mass in that snapshot ra m ax. Masses are in units of 
10 10 M o . 



merge into larger haloes. Note that Kly pin et al.| ( |2003] ) conclude 
the opposite. This is because they normalise the power spectra of 
the different cosmologies at z = instead of at z = 1089 as we 
do. In we discuss this matter further. 

The other panels show (clockwise, starting at the top right) the 
mean halo mass (m), the maximum halo mass in the sample m max 
and the standard deviation in the mass distribution cr m . These pan- 
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els show that at every redshift the haloes in the SUGRA cosmology 
are less massive than those of the others, implying that evolution 
of structure in the SUGRA universe is slowest. See De Boni et alTl 
(201 1) for a thorough analysis of the halo properties and their pos- 
sible direct use as probes of the nature of dark energy. 



4 VOID IDENTIFICATION & ANALYSIS 

Our void analysis of the particle and halo samples involves a few 
steps. In this section we will describe these in some detail, along 
with a schedule of the analysis setup. 

The fundamental step of the entire analysis is that of finding 
and defining voids. For this we use the Watershed Void Finder 
(Plat en, van de Weygaert & J ones 2007). Having identified the 
voids, we determine the shapes and sizes of the voids. Subse- 
quently, we need to identify and evaluate the possible systematic 
effects that play a role in the inferred void properties. Also, for the 
interpretation of the significance of the results, we have to take into 
account the cosmic variance. 



4.1 Watershed Void Finder 

Depending on your definition, voids can make up from 13 to 100% 
of the total volume. Evidently, there is no unanimous agreement 
on the definition of void. It is one of the reasons why there are 
a wide range of different void finding algorithms, often yielding 
vastly different results (see e.g. |Colberg et al.|2008| >. For our sta- 
tistical analysis we use the void finder developed by our group, 
the Watershed Void Finder (WVF), introduced by Plat en, van de| 
Wey gaert & Jones | ( |2007| >. It is a largely parameter free formal- 
ism that manages to outline a void region, independently of its size 
and shape, and thus ideally suited for an objective statistical study. 
Moreover, the study by |Platen, van de Weygaert & Jon es (2008) 
showed that the voids identified by WVF are indeed closely related 
to the physical structure of the mass distribution: their orientation 
is intimately coupled to the tidal field, closely following the related 
theoretical predictions. 

The Watershed Void Finder (WVF) is an implementation of 
the Watershed Transform for segmentation of images of the galaxy 
and matter distribution into distinct regions and objects and the sub- 
sequent identification of voids. The basic idea behind the watershed 
transform finds its origin in geophysics. It delineates the bound- 
aries of the separate domains, the basins, into which yields of e.g. 
rainfall will collect. The analogy with the cosmological context is 
straightforward: voids are to be identified with the basins, while the 
filaments and walls of the cosmic web are the ridges separating the 
voids from each other. 

With respect to the other void finders the watershed algorithm 
has several advantages. Because it identifies a void segment on the 
basis of the crests in a density field surrounding a density minimum 
it is able to trace the void boundary even though it has a distorted 
and twisted shape. Also, because the contours around well chosen 
minima are by definition closed, the transform is not sensitive to 
local protrusions between two adjacent voids. The main advantage 
of the WVF is that for an ideally smoothed density field, it is able 
to find voids in an entirely parameter free fashion. 

The WVF consists of eight steps, which are extensively out- 
lined in Platen, van de Weygaert & Jones (2007). For the success 
of the WVF it is of utmost importance that the density field re- 
tains its morphological character. To this end, the two essential first 
steps relate directly to the Delaunay Tessellation Field Estimator 




(DTFE), which guarantees the correct representation of the hierar- 
chical nature, the weblike morphology dominated by filaments and 
walls, and the presence of voids ( Schaap & van de Weygaert \ 
|van d e Weygaert & Platen 201 1 Cautun & van de Wey gaertp 
The obtained densities are interpolated onto a regular grid of 384^ 
grid cells. 

Because in and around low-density void regions the raw den- 
sity field is characterised by a considerable level of noise, a second 
essential step suppresses the noise by means of a properly defined 
filtering operation. In its original definition, we invoked an adaptive 
smoothing algorithm, median filtering, which in a consecutive se- 
quence of repetitive steps determines the median of densities within 
the contiguous Voronoi cell surrounding a point. The determination 
of the median density of the natural neighbours turns out to define 
a stable and asymptotically converging smooth density field fit for 
a proper watershed segmentation. An alternative filtering operation 
is one in which we smooth the field by Gaussian filters of a fixed 
scale. The WVF segmentation of the filtered field produces the void 
population corresponding to that one particular scale. In principle, 
one can combine this with the WVF void segmentations at other 
scales. By carefully examining and linking small scale void struc- 
ture and the large scale voids one may infer the multi- scale structure 
of the hierarchically evolved void population ( [Sheth & van de Wey^| 
|gaert|2004l|Aragon-Calvo & Szalay|2012| ). In our study, we have 
opted for this last filtering operation, as it will allow us to explicitly 
assess the scale dependence of the void shape measurements. This 
will enable us to find explicit scale criteria for a successful dark 
energy exploration based on the void population. 

The subsequent central step of the WVF formalism consists 
of the application of the discrete watershed transform on this adap- 
tively filtered density field. As a result of noise in the (discretely) 
sampled density field, the WVF segmentation includes artefacts. To 
a first order these can be classified as false splits and false mergers 
of void-patches. To deal with this these effects, we include a step 
of removal of segmentation boundaries whose density is lower than 
some user-defined density threshold. 

To obtain an impression of the application of the WVF for- 
malism, for two different scales figure [6] illustrates the resulting 
watershed segmentation of the density field in the ACDM N-body 
simulation. In the four bottom panels, the top two panels have 
been smoothed with a Gaussian filter radius Rf — 1.5 /i _1 Mpc, 
the bottom two panels with Rf = 6.0 h~ 1 Mpc. In both cases 
the void boundaries follow the filamentary structures quite accu- 
rately and the irregular shapes of the voids are well preserved. 
Also, comparison between the watershed void boundaries of the 
R f = 1.5 h^Mpc and R f = 6.0 /i _1 Mpc filtered fields shows 
the loss of void substructure on the larger scale. 

A related tessellation-based method for void identification, 
ZOB OV (Neyrinck|2 008), does yield similar results as WVF ( [CoT] 
|berg et al" |2008| ). It demonstrates the successful application of 
tessellation-based techniques to identify structures within the cos- 
mic matter distribution. 



4.2 Void shapes 

After determining the void distribution, we calculate properties of 
each of the voids. This includes some fits of shape parameters. In 
what follows we represent the density field on uniform grid cov- 
ering the volume of interest. A void consists of the group of grid 
cells that WVF has identified as belonging to a particular watershed 
basin. The geometric centre of a void is defined as the (volume) av- 
erage of the void's grid cell, or voxel, positions. The volume of the 
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Figure 6. Top: a density field slice of thickness 0.25 h~ 1 Mpc from the ACDM simulation. This slice of 300 by 150 h~ 1 Mpc is the same as those used 
in the bottom panels. Colour represents density as in e.g. figures [7] and [9] Bottom left two: the corresponding distribution of voids (void borders in black 
lines, the different colours do not represent anything physical). For the top panel we used a Gaussian filter radius Rf = 1.5 h~ 1 Mpc, for the bottom panel 
Rf = 6.0 h~ 1 Mpc and the 1.5 h~ 1 Mpc ones are transparently inset. Bottom right two: a random selection of ellipsoid fits (yellow) overlaid on the density 
field (now in grayscale), again at two radii 1.5 and 6.0 h~ 1 Mpc. 



voids is determined by simply adding the volumes of the voxels 
that define the void. 

To determine the shape of a void, we fit its volume by an el- 
lipsoid. We assume the density of the void to be uniform. The ap- 
proximation of voids as a homogeneous ellipsoidal region is a more 
than adequate first order approximation for the interior of voids in 
a wide range of cosmologies, justified by our understanding of the 
formation of voids around minima in the primordial density field 
(e.g. |Icke|| 1984[ |van de Weygaert & van Kampen][T993l |Sheth &| 
van de Weygae rt|2004[ |Shandarin et al. 2006}. To some extent, it 
is a considerably better description for voids than it is for over- 
dense regions. Overdense regions contract into more compact and, 



hence, steeper density peaks, so that the area in which the ellip- 
soidal model represents a reasonable approximation will continu- 
ously shrink. By contrast, the region where the approximation by a 
homogeneous ellipsoid is valid grows along with the void's expan- 
sion. While voids expand, their interior gets drained of matter and 
develops a flat "bucket- s haped" density profile | van de Weygaert &| 
|van Kampe"n[ ( |1993| ); |Sheth & van de Weygaert| ( |2004| ). Hence, the 
void's natural tendency is to evolve into expanding regions of a 
nearly uniform density. The approximation is restricted to the inte- 
rior and fails at the void's outer fringes because of its neglect of the 
domineering role of surrounding material, such as the sweeping up 
of matter and the encounter with neighbouring features. 
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The homogeneous ellipsoidal shape allows us to focus en- 
tirely on the geometrical properties of the void and avoid possible 
complications introduced by the overdense regions in and around 
the void. In practice, we proceed as follows. We first calculate the 
void's inertia tensor Uj\ 

lij = _ X kiXkj) , (11) 

k 

where we sum over all cells k belonging to the void. In the above 
expression, Xk = (xfci, Xfc2, Xkz) is the distance vector of the /c-th 
void cell to the void's centre, and Sij is the Kronecker delta. 

The ellipsoidal fit to the void is taken to be the one that would 
have the same inertia tensor hj as the void. 

The size, shape and orientation of the ellipsoid are inferred 
from the eigenvalues and the eigenvectors of the inertia tensor Uj . 
The relevant properties in this work are the ellipsoid's semi-axes a, 
b and c (a ^ b ^ c), which are connected to the eigenvalues of the 
inertia tensor as: 

a 2 = \(h + h-h) , (12) 
b 2 = |(/ 3 + /i-/ 2 ) , (13) 
c 2 = ^{h+h-h) • (14) 

Subsequently, we can characterise the shape of the void by three 
parameters: the ellipticit)[^] e or conversely the sphericity s, the 
oblateness (flattening) p and the prolateness q. These quantities are 
defined as: 

1 C ° h A C n<\ 

e—1 , s — -. V — and q — t ■ (15) 

a a a b 

The ellipsoid's volume is V = ^Trr 3 , where r = Vabc is the void's 
effective radius. Note that this volume is slightly different from the 
real void volume due to irregularity of the void's borders. 

To obtain an idea of the resulting void ellipsoids, the bottom 
righthand panels of figure[6]depict a random selection of void ellip- 
soids, superimposed on the density field. The top frame shows the 
ellipsoids corresponding to the WVF voids obtained after smooth- 
ing the density field by 1.5 /i _1 Mpc, while the bottom frame shows 
the void ellipsoids corresponding to the 6 h~ 1 Mpc smoothed den- 
sity field. In general, there is a good correspondence between the 
size and shape of the ellipsoids and the voids in the underlying den- 
sity field. Evidently, the fits are rarely perfect as the void boundaries 
tend to have a rather irregular shape. One consequence of this is that 
the ellipse volumes differ slightly, by an average factor of ~ 1.08, 
from the actual void volumes, determined by adding the volumes 
of the grid cells that WVF identified as belonging to its interior. 

4.3 Systematics 

The measurement of void shapes in the observational reality, and 
the corresponding statistical analysis, is subject to a number of arte- 
facts and systematic effects. The three main effects that we include 
in our analysis are galaxy biasing, redshift distortions, and cosmic 
variance. 



1 For consistency, we follow |Park & Lee (2007) and others in using the 
term ellipticity for this quantity. This term may not be the most adequate, 
as its semantics imply a general description of the shape of an ellipsoid, 
whereas the total shape of an ellipsoid needs at least two parameters. "As- 
phericity" may therefore have been a better term. 



4. 3. 1 Discrete and diluted samples 

By far the most complex effect is that of galaxy biasing. Within the 
context of this study, we concentrate on the consequences of the 
discrete and diluted nature of galaxy samples with respect to the 
underlying mass distribution, and the biasing of the halo popula- 
tions. 

We will not address the effects of the far more complex is- 
sue of biasing on the basis of intrinsic galaxy properties, and re- 
fer to an upcoming study in which we systematically address the 
properties of voids in galaxy populations modelled by different 
semi- analytical models of galaxy formation. Here we address the 
effects of the discrete and diluted nature of the galaxy distribution 
by means of random subsets of simulation particles of the N-body 
simulations ( §3.1| >, in combination with samples of haloes from the 
same simulations. 

4.3.2 Redshift distortions 

The main distorting effect that we consider is that resulting from 
redshift distortions due to the peculiar velocity of matter, haloes 
and/or galaxies. Given that galaxy redshift survey data distances are 
purely based on redshifts, which inevitably introduces an error due 
to the considerable contribution of radial velocities to the redshift. 

As a result, the shapes of voids are expected to get system- 
atically distorted. The void will apparently expand and, hence, in 
redshift space they will appear more elongated along the radial di- 
rection. This will result in a redshift space void distribution that is 
expected to be systematically shifted to slightly higher values of 
the ellipticity e. Also, the orientation of the voids in redshift space 
will be somewhat anisotropic, with a slight tendency to be directed 
along the radial direction. The effect is rather small, and only be- 
comes apparent when considering a large and statistically represen- 
tative sample of voids (for a thorough recent analytical treatment, 
see |Shoji & Lee|2012| >. 

In this study we restrict ourselves to estimating the (cosmic) 
variance in ellipticity due to the fact that we do not know the true 
direction and magnitude of the peculiar velocity of dark matter 
and haloes. We will take this into account as an extra error term 
when assessing the significance of systematic differences between 
the void shape distributions in different cosmologies. To get these 
estimates, for each of the five cosmologies we generate eight red- 
shift space realisations of the matter and halo distributions. Each 
of the redshift space realisations is defined with respect to a ran- 
domly chosen location in the simulation box, the "observer" lo- 
cation. Different observers see different redshift spaces due to the 
differences in radial velocities with respect to the observer. In each 
of the redshift space realisations we identify the void population 
and determine the void shape distribution using the same tools as 
in the underlying real space realisation. 

Following the choice of an observer's position, the redshift 
space realisation is computed from a given simulation by trans- 
forming the particle locations to redshift space with respect to the 
"observer". Locating the observer at the centre of the box (and 
translating the particle and halo positions accordingly), the mea- 
sured redshift of the particle (or halo) is 

Z = (1 + 3 H )(1 + Z r ) ~ 1 ~ ZH + Z r ~ — + — , (16) 

c c 

where v r is the radial component of its peculiar velocity with re- 
spect to the observer, 

v r = v pec • r . (17) 
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Figure 7. Redshift space comparison. Shown are "redshift space" (see ^4.3.2) density slices of thickness 0.25 h 1 Mpc of the simulation boxes of the ACDM 
(left) and SUGRA (right) models at z = 0. The observer is placed at the centre. Density values are related to the colours by the bar on the right. 



In this expression, zh is the cosmological redshift of the particle or 
halo due to the Hubble expansion. Its resulting redshift space posi- 
tion is then given by (r z _ spa ce, 0, </>), with and (j) its sky position 
as seen from the observer and 

r z - S pace = ZC/H . (18) 

Examples of resulting redshift space realisations are shown in 
figure [7] They show the redshift space structure for a ACDM and 
a SUGRA cosmology, for the same observer's position and time. 
Large clusters are transformed into the well known Fingers of God 
(Jackson 1972), whereas voids in redshift space appear to be some- 
what larger and emptier than in the corresponding real space dis- 
tribution. Differences between the two models are similar to those 
seen in real space. The most prominent difference is that of the 
resulting Fingers of God. In the SUGRA cosmology these clus- 
ters appear to be significantly shorter than their equivalents in the 
ACDM cosmology. This indicates a lower internal velocity disper- 
sion of clusters in the SUGRA cosmology, a manifestation of the 
slower evolution and lower mass of clusters (see §5.1| >. 



4.3.3 Cosmic variance 

We additionally include the error introduced by cosmic variance. 
We calculate the amount of variation in the evolution of the main 
void ellipticity by taking eight different random phase realisations 
of our low resolution simulations and calculating the standard devi- 
ations in the eight obtained mean ellipticity values at all redshifts. 

The value of the deviation is quadratically added to the esti- 
mated redshift distortion. These combine to the square of the total 
estimated "error", which we use in the interpretation of the results 
presented in the figures of ^6] 



4.4 Analysis setup 

In our analysis, we compare the void shape results from the dif- 
ferent cosmological samples to each other and to previous results 
from the literature. The analysis involves the following steps: 

• We first compare our results to those found in the literature. 
The high resolution (768 3 ) simulations are used to see whether 
the mean void ellipticity indeed evolves as predicted, and increases 
with decreasing redshift. 

• On the basis of the low resolution (256 3 ) data, we assess 
whether we can distinguish between the dark energy models by 
comparing their redshift evolution of the mean void ellipticity, 
(e)(z). 

• We compare the low and high resolution results to check for 
effects of mass resolution. There is no physical reason for there to 
be a resolution effect, so any effect will be due to the methods used. 

• For each of the simulations, we compare the results of the av- 
erage void shapes obtained for the pure DM particle distributions 
with those obtained from the corresponding DM halo populations. 
To a large extent the DM halo distribution reflects the galaxy dis- 
tribution. In terms of estimating the effect of its discreteness and 
sparseness, the DM haloes are largely representative. Results on 
DM haloes may therefore be used to answer the question of whether 
or not the inferred void shape evolution represents a significantly 
discriminative probe of dark energy in the case of real galaxy red- 
shift surveys. 

We find marked differences in (e)(z) between DM haloes and DM 
particles. To identify the origin of these differences we run the fol- 
lowing additional tests: 

• We compare (e)(z) for the DM haloes and samples of the 
same number of randomly sampled simulation particles from the 
high resolution simulations. In case the results of these two differ- 
ent samples are the same, we may conclude that the void geometry 
of the halo distribution does not intrinsically differ from that of 
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Figure 8. Void evolution in ACDM. Part of the density slices (range: x ~ [47, 203], y ~ [47, 188] /i -1 Mpc) of thickness 0.25 /i _1 Mpc of the simulation 
boxes of the ACDM model at redshifts (from top to bottom) z = 2.0, 1.0, 0.5, 0.25, 0.1 and z = are shown. The colors were chosen to highlight the 
evolution of voids, especially in terms of their substructure: dark blue indicates the lowest density regions, green and white indicate the higher density regions 
(see colour bar on the right). 
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Figure 9. Evolution of structure in different cosmologies. Part of the density slices (range: x ~ [47, 203], y ~ [47, 188] h~ 1 Mpc) of thickness 0.25 h~ 1 Mpc 
of the simulation boxes of the ACDM, SUGRA, RP, EQn and EQp models (from top to bottom) at z = are shown. Cluster regions are yellow (p/po > 200) 
and voids are black (p/po < 0.2), as also indicated in the colour bar to the right. The zoom-in boxes (top right) focus on comparable cluster regions in the 
ACDM and SUGRA simulation. 
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the dark matter distribution. It would imply that the differences be- 
tween the void shapes found in the haloes and in the dark matter 
(particle) distribution are mainly the result of the sparsity of the 
halo population. On the other hand, if there are notable differences 
between the halo voids and the dark matter particle voids, the bi- 
asing of haloes - and by implication also those of galaxies - is an 
important factor. 

• As an additional - strong - test of possible systematic differ- 
ences between the spatial distribution of haloes and dark matter 
particles, we compare the random samples of dark matter particles 
to the unweighted halo distribution. The latter consists of the same 
haloes as in the regular halo sample. However, instead of includ- 
ing their true masses, we assign the same mass to each halo. If 
the mean void size (e) (z) differs even more in this test than in the 
previous one, the void population in the halo distribution has to be 
significantly different from that in the dark distribution. This would 
imply the inescapable conclusion that biasing is an important cause 
for the differences in (e) (z). 



5 RESULTS: VOID POPULATION CHARACTERISTICS 

In this section we will first present some qualitative results on the 
evolution of the WVF identified void population and discuss some 
overall quantitative characteristics of the void populations. 



5.1 Voids in the dark matter field: visual impression 

The evolving dark matter distribution in the ACDM cosmology is 
shown in figure [8] We can see the evolution of several voids (dark 
blueish regions). We roughly define these as regions with a density 
contrast of p/po < 0.2. 

Overall, we observe that the most prominent evolution of the 
void population takes place between redshifts 2 and 0.5. This agrees 
with |Huterer & Turn er (2001 ), who argued that the redshift range 



0-2 ^ z < 2 is the most promising for probing w(z) (see also Hell 



wing, Juszkiewicz & van de Weygaert 2010). The voids quickly 
grow from initial seeds with radii of about 1 /i _1 Mpc, to voids 
ranging in size from 2 to 30 /i _1 Mpc. We also see that in this pro- 
cess the number of voids decreases due to the merging of voids, a 
manifestation of the hierarchical evolution of voids (She th & van de| 
|Weygaert|2004l|Aragon-Calvo & Szalay|2012) . At the later time- 
steps, we may still discern a lot of substructure in voids. These are 
the remnants of the same structural evolution: even billions of years 
after the voids merge, their outline remains visible as tenuous un- 
derdense features in the realm of the emerging void (see S heth &| 
|van d e Weygaert 2004 ). This is clearly visible in the bottom pan- 
els, corresponding to z — 0.1 and z = 0.0, where we may clearly 
discern walls in the larger voids. 

In figure [2] we compared the evolution of the void population 
in the ACDM simulations with that in the other four cosmological 
models. The most outstanding differences are those between the 
ACDM cosmology on the one hand, and the SUGRA cosmology at 
the other extreme. 

The large scale structure in the SUGRA cosmology at z — 
is less evolved than that in the ACDM cosmology. One of the man- 
ifestations of this concerns the voids, whose size is smaller. Also, 
the filaments in the SUGRA cosmology are more diffuse, and clus- 
ters are less pronounced and clumpier. The lower level of evolution 
in the SUGRA model ties in with the higher value of w over a large 
fraction of its cosmic evolution. The closer the value of the dark 
energy equation of state w is to its upper limit of —1/3, the shorter 



the cosmic timescale over which structure can evolve. The SUGRA 
universe, for instance, is effectively ~ 560 Myr younger at z = 
than the ACDM universe. One can also observe this from the mea- 
sured values of as- As shown in figure [3] at each redshift SUGRA 
has a lower value of as than that found in the other cosmologies, 
while ACDM has the largest amplitude of mass fluctuations and 
structure. 

The RP, EQn and EQp dark energy models differ less dramati- 
cally from the ACDM model. All largely follow the same evolution, 
and it is more difficult to visually distinguish the resulting struc- 
tures from those seen in the ACDM cosmology than it was for the 
SUGRA cosmology. The void regions observed in the correspond- 
ing frames of figure|9]have a large degree of similarity in shape and 
size to those found in the ACDM simulation. 

To emphasize the large difference between the ACDM cos- 
mology on the one hand, and the SUGRA cosmology on the other 
hand, we zoom in on a central cluster region (see zoom-ins in fig- 
ure [9]). When comparing these objects, we should note that the 
corresponding simulations started from the same initial conditions, 
so that any difference in structure, morphology and dynamics is 
a direct reflection of the influence of dark energy on the structure 
formation process. Not only do we find that the ACDM cluster is 
more compact and centrally concentrated, but also that the mor- 
phology of the surrounding matter distribution is substantially dif- 
ferent. The weblike filamentary structures in the SUGRA cosmol- 
ogy find themselves in a dynamically younger state. The number of 
filaments is higher, while they are thinner and less well defined. In 
the ACDM situation, we find the central cluster embedded in a web 
of a few strong and well defined filaments. 



5.2 Void sizes and shapes 

Using the void population obtained by means of the WVF, the size 
and shape of each of the voids is calculated following the descrip- 
tion in 



5.2.1 Different cosmologies at z = 

In figure [TO] we plot the (effective) void radii r and ellipticities e in 
the five cosmologies, at z = 0. These are all computed from the 
high resolution simulations. The functions shown in the plots are 
the normalised probability density functions, for which 



f 

J —c 



f(r)dr = 1 ; 



f 

J —c 



f(e)de = 1 . 



(19) 



These results were obtained at a smoothing radius of Rf = 
1.5 /i _1 Mpc. 

One conclusion from this is that the effects we are trying to 
measure are rather subtle. To uncover systematic trends, we there- 
fore will have to direct our attention to moments of the ellipticity 
distribution and average over the entire sample. In this respect, we 
follow earlier studies in concentrating on the mean ellipticity. 



5.2.2 Scale dependence 

One of the determining factors in our analysis is the scale at which 
we study the void population. 

The effect of the smoothing radius on the size distribution of 
the voids is shown in figure |ll(a)| We see that the size of the de- 
tected voids increases systematically with smoothing scale. This is 
because the small scale void boundaries are smoothed out, resulting 
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Figure 10. Void radius and ellipticity PDFs (R f = 1.5 h~ 1 Mpc) in the 
different dark energy model simulations at z = 0. 
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Figure 11. Void property PDFs for different WVF smoothing scales Rf. 



in a sample of larger void regions (see figure[6]). The distribution is 
peaked around a characteristic void size. The void population is 
marked by a large number of voids around this characteristic void 
size, with hardly any smaller voids and a strongly declining number 
of larger voids. This distribution agrees entirely with that predicted 
within the context of the excursion- set description of the hierarchi- 
cal void evolution ( Sheth & van de Weygaert 2004). As is born out 
by the figure and as expected, the peak size - as well as the en- 
tire distribution - systematically shift towards higher values as the 
smoothing scale increases. Indications towards this behaviour have 
independently also been found in other studies. One example is the 
recent study by |Einasto et al.| ( |2011| ), who used a different void 
finder and different void-delineating objects. They noted a system- 
atic shift towards higher (mean) void radii with increasing lower 
threshold mass of the defining objects. This interesting parallel was 
to be expected, as the high mass objects will form preferentially in 
the high density regions that will dominate the cosmic web on the 
largest (smoothing) scales. 

Of more immediate interest for our program is the behaviour 
of void shapes as a function of scale. The void shape distribution 
turns out to vary only mildly with the scale of the density field, 
as shown in figure [TT(b)] To a first approximation, the shape dis- 
tribution of the void population appears to have a nearly scale in- 
variant character. There is a slight tendency for the voids on the 
largest scales of 6.0 and 8.0 h~ 1 Mpc to be slightly more spheri- 
cal. However, to a large extent this may be ascribed to the fact that it 
concerns samples with a considerably lower number of objects and 
hence with a larger statistical uncertainty. Nonetheless, the nearly 



scale independent behaviour of the void shape distribution is an in- 
teresting and highly relevant result. It implies that - in the case of 
the dark matter voids - we can restrict ourselves to an evaluation 
of the shapes of voids on a scale of Rf — 1.5 /i _1 Mpc, as this 
would contain all necessary shape information and guarantee a sta- 
tistically optimal result. 

5.2.3 Overdense void boundaries 

One feature of the WVF algorithm that requires some attention is 
the fact that, in its pure form, it treats the entire interior of each 
void basin as belonging to a void. The only regions which are not 
included are the boundary grid cells where two voids meet. We in- 
vestigated whether the inclusion of overdense regions at the bound- 
aries of detected voids might cause artefacts in the inferred void 
shape distribution. Only the inner parts of the void are properly de- 
scribed by a homogeneous ellipsoidal model, [^] and are therefore 
expected to adhere better to the ellipsoidal fit of our analysis. 

By removing the overdense boundary regions, about 14% of 
the volume in voids is eliminated. The shape analysis based on 
these voids does indeed yield slightly different results. One finding 
is that the mean ellipticities tend to be slightly higher. However, 
the offset is quite comparable for all simulations, at all redshifts. 
There are no changes with respect to any of the qualitative results 

2 Because the "role of surrounding material [will dominate], through the 
sweeping up of matter and the encounter with neighbouring features" ( |van| 
|de Weygaert & Platen|20TT) . 
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Figure 12. Void ellipticity PDFs (R f = 1.5 /i _1 Mpc) in the ACDM 
model in real space and in 8 different redshift spaces. 



discussed in the following sections. We therefore decided to keep 
to the basic WVF algorithm. 



5.3 Voids in redshift space 

In figure [T2| we compare the ellipticities of voids in the ACDM 
dark matter distribution with that in eight redshift space realisations 
within the same high resolution simulation. We find a rather consis- 
tent offset of the ellipticity distribution of around Ae ^ 0.014. The 
redshift space voids have, as expected, a systematically more elon- 
gated shape. In the other cosmologies, the situation is comparable, 
but the effect is somewhat smaller. 

In addition to the systematic shift, which represents the sys- 
tematic redshift space distortion, there is also a sizeable scatter 
between the different redshift space realisations. This is the error 
which is introduced by local variations in the contribution of pe- 
culiar velocities to the redshift of objects. It behaves like an extra 
contribution to the cosmic variance. The scatter is also a function 
of filter scale Rf, varying from 0.00076 in the case of voids on a 
scale of Rf = 1.5 /i -1 Mpc to 0.0012 for voids on a larger scale 
of Rf — 3.0 /i _1 Mpc. In the subsequent analysis we restrict our- 
selves to voids on these scales. 
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Figure 13. Void radius and ellipticity PDFs of voids in the dark matter and 
halo distributions, at different scales. 



Both their size and shape distribution are, within acceptable limits, 
similar. 



5.4 Voids in the halo distribution 

When comparing the spatial halo distribution in figure [4] with the 
underlying dark matter distribution, we observed the substantial 
loss in spatial resolution. The halo distribution, even at z = 0, is 
much sparser than the dark matter particle distributions in the sim- 
ulations. The diluted halo sample may therefore not be expected 
to trace the fine features visible in the matter distribution, and cer- 
tainly will not be able to accurately sample the small scale void 
population. 

In how far this affects our study may be appreciated from fig- 
13 At a scale of Rf = 1.5 /i -1 Mpc, the haloes completely 
fail to reproduce the size distribution of the voids in the underly- 
ing void distribution (figure [l3(a)| ). This may be hardly surprising: 
haloes cannot resolve these structural features. Along with this, we 
also see that the halo distribution must be used with care when in- 
vestigating the void ellipticity distribution (figure [13(b)] ). 

On the other hand, at a large scale, like Rf = 8 /i _1 Mpc, 
we find complete accordance between the void population in the 
dark matter distribution and the voids found in the halo distribution. 



5.4.1 Scale dependence and sample size 

In the generic observational situation, where the voids are traced by 
a discrete and dilute galaxy distribution, the scale dependence of 
the void analysis involves two opposing effects. On the one hand, 
it is easier to reliably trace the outline and measure the shape of 
larger scale voids. On the other hand, the number of voids traced 
by the object sample within a given volume of observed space will 
strongly decrease as we consider voids on larger scales. The errors 
in the measured statistical moments will increase accordingly. 

Hence, it is necessary to find a compromise between both re- 
quirements. However, the sampling density of objects may get so 
sparse that it becomes unfeasible to trace enough voids on a suffi- 
ciently large scale. In that situation, it will be impossible to use the 
void population for an attempt to measure the dark energy equation 
of state as any such measurement will suffer from inherently large 
errors. 

One possible improvement of this will be to substantially in- 
crease the survey volume. This would enable to obtain more re- 
liable and representative statistics for the larger number of small 
scale voids. 
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Figure 14. Void ellipticity PDFs (R f = 1.5 /i _1 Mpc) in the ACDM 
model simulations at different redshifts. The bottom panel shows the ratio 
of the PDFs at z ^ to that at z = 0. 

6 RESULTS: VOID SHAPE EVOLUTION 

Having established the basic size and shape properties of the void 
population in the different simulated cosmologies, we arrive at the 
examination of the time evolution of the mean ellipticity, 

(e)(z) = J deef(e,z) (20) 

and its relation with the character of the dark energy in the corre- 
sponding cosmology. 

6.1 Basic results vs literature 

The void ellipticity distribution in the ACDM cosmology evolves 
systematically with redshift. This is clearly visible in figure [14] 
where we find a gradual shift of the ellipticity distribution towards 
higher values of e as time proceeds. In other words, the ellipticity 
of voids is expected to decrease towards higher redshifts. This is 
entirely in line with the expected generic behaviour, as described 
by e.g. |Lee&Park| ( [2009t . 

The evolutionary trend of the mean ellipticity (e) follows the 
general trend of |Lee &~ Park (2009 ), as is clearly shown in fig- 
ure |15(a)| The mean ellipticity of the voids increases with time, 
with voids being less elongated towards higher redshifts. In one re- 
spect our results differ with those obtained by ( |Lee & Park 2009 ; 
|Lavau x & Wandelt 2010 ): the voids in our simulations do not have 
equal ellipticities at z — 0. This is a result of the difference in 
normalisation between our simulations. By normalising our simu- 
lations, via as, at the recombination redshift z — 1089, the level 
of structure formation at z — between the different cosmologies 
may be expected differ. The differences in mean void shapes is one 
particular manifestation. 

When comparing the ACDM and SUGRA results for the el- 



lipticity evolution, we find that there is a significant redshift range 
over which we can clearly distinguish between the void ellipticities 
in the different cosmologies. 

6.2 Simulation resolution 

When comparing the inferred mean void ellipticity evolution in the 
low resolution 25 6 3 ACDM simulations and the high resolution 
768 3 ACDM simulations, we find that there is a slight resolution 
effect. We analysed the two simulations, which have exactly the 
same initial conditions - in terms of mode amplitudes and phases 
- and cosmological parameters. Figure [T5(a)| reveals the difference 
between the high resolution (yellow) and low resolution (blue) sim- 
ulation. The only explanation is the difference in mass resolution of 
the simulations. 

To confirm that the difference is purely a resolution effect, 
we compared the low and high resolution simulations at large fil- 
ter radii Rf, to smooth out the small scale differences caused by 
the different resolutions. Indeed, we see in figure |15(e)| that at large 
scales, the results for the high and low resolution simulations con- 
verge. 

We also find that at these large scales the power of the analysis 
to discriminate between different dark energy cosmologies is lost. 
The figure reveals that it is not possible to distinguish significantly 
between the ellipticity evolution curves obtained from the SUGRA 
simulation and the corresponding ACDM curves. This strengthens 
our choice to opt for a study of the void population at filter scales 
of R f = l.bh^MpcmdRf = 3.0 /i _1 Mpc. 

6.3 Voids in the halo distribution 

Given that the resolution of the DM halo sample, of around 
560,000-580,000 particles, is considerably lower than that in the 
low resolution DM particle sample we may not be surprised to find 
that it is difficult to find any significant evolutionary trends in the 
corresponding void population. 

We indeed find that it is not possible to detect a decrease in 
mean ellipticity of the void population. Rather, we find that the halo 
void ellipticities stay remarkably constant up to at least z = 1 and 
rise steeply towards higher redshifts (see figure [T5(b)] ). The latter 
is probably a result of the (strongly) decreasing number of haloes 
at higher redshifts and the corresponding poorer sampling of the 
underlying full density field. Any hope seems lost to discriminate 
between different dark energy models on the basis of the measured 
void shape parameters. 

6.4 Sparsity effect 

A key question for understanding the inability of haloes to repro- 
duce the ellipticity evolution of voids is whether this is mostly 
an effect of the discreteness and sparseness of the halo sample or 
whether intrinsic biasing is also at play. 

To this end, we take a random subsample of particles from the 
dark matter particle distributions and repeat the analysis. The num- 
ber of randomly sampled particles is taken to be equal to the halo 
population. Bootstrapping errors are used to obtain error estimates 
of the inferred mean void shapes. They are the standard deviations 
of the eight mean ellipticities that were obtained from eight differ- 
ent random subsets of the high resolution DM particle set. 

Figure [15(c)] shows that the (lack of) void shape evolution in 
the diluted random DM particle distribution largely agrees with the 
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(a) Mean ellipticity as a function of redshift. ACDM high and low resolution 
simulations and the SUGRA low resolution one. 
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(b) Mean ellipticity as a function of redshift. ACDM and SUGRA low reso- 
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(e) Mean ellipticity as a function of redshift. The low and high resolution 
ACDM results are shown at higher filter radii. The low and high res simula- 
tions converge at these higher filter radii, because the small scale differences 
are smoothed out. 

Figure 15. Mean ellipticity as a function of redshift, (e) (z). In the four plots above, in the right part the z-axis is in logarithmic scale and the left part is in 
linear scale (so we can show z = as well). 



© 2012 RAS, MNRAS 000,[l]{22] 



The darkness that shaped the void 19 



results of the halo void study. This implies that the deviation of the 
halo void shapes from the shapes of voids in the high resolution 
dark matter simulations is to be ascribed to the sparsity of the halo 
population. 

However, we also observe some additional differences. On a 
scale of Rf = 3 h~ 1 Mpc, the void shape evolution in the ran- 
dom subsample appears to differ significantly from that of the halo 
void ellipticity curves. Moreover, the increase of ellipticity with 
redshift that is observed in the halo void sample (at a scale of 
Rf = 1.5 /i -1 Mpc), is not reproduced by the voids in the subsam- 
pled particle distribution (at Rf = 3 /i _1 Mpc). This argues for the 
influence of additional effects, in particular that of the spatial bias 
of haloes. 

This conclusion is confirmed when considering the results ob- 
tained for the void population in the unweighted halo distribution 
(see §4.4| ). The differences with the random subsample voids be- 
come slightly larger than in the case of the reg ular we ighted halo 
voids, especially at Rf = 3 h~ 1 Mpc (see figure 15(d)| . 

The inescapable conclusion appears to be that the spatial bias 
of the halo population, and all accompanying complications, is in- 
deed a factor of significant importance when seeking to use the void 
population as tracer of the dark energy equation of state. 



7 SHAPES AND CLUSTERING: THE a 8 DEGENERACY 

We looked into the possibility that the differences between the 
model ellipticities are not primarily caused by differences in the 
equation of state parameter w(z) of dark energy, but by some other 
cosmological parameter. We have come to the interesting conclu- 
sion that the one exclusive and dominant factor is that of the level 
of clustering and structure development, expressed by as. In other 
words, the void ellipticity distribution is a manifestation of the level 
of clustering of matter. 

At the present epoch, the different dark energy models have 
different values of as. This is the result of the different structure 
growth rates between the different cosmologies, all starting from 
the same primordial density field whose amplitude has been nor- 
malised at the time of last scattering. The value of the amplitude is 
fixed by the value of as measured from WMAP data. 

In figure [16] we see that there is a strong correlation between 
as and the mean ellipticity, independent of dark energy model or 
redshift. Both lines consist of as measurements at redshifts of 0, 
0.1, 0.25, 0.51, 1.0 and 2.04. If they had been redshift dependent, 
we should have been able to distinguish e.g. the point at the middle 
of the SUGRA line (which corresponds to z ~ 0.5) from the point 
at the same location of the ACDM line (which corresponds to z ~ 
1). The strong correlation between as and ellipticity is especially 
surprising when we realise that the ACDM and SUGRA models 
differ substantially in almost ever other aspect that we investigated. 
Yet, there is little to no difference between these two models in 
figure [16] 

It seems then that the main cause of the differences between 
these models is, in fact, as. Of course, in this case, the difference 
in dark energy models, i.e. in w(z), is the underlying cause of the 
differences in as. w(z) exerts its influence on as through D(z) 
(equation [9}, which in turn contains a factor H(z), thus realising 
the coupling to w(z) (equation [5]). 
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Figure 16. Mean ellipticities versus as. The different lines show the low 
resolution simulations of the five different cosmologies. The lines consist 
of cr 8 values at redshifts 0, 0.1, 0.25, 0.51, 1.0 and 2.04. 



8 CONCLUSIONS AND DISCUSSION 

We have investigated the suggestion that the evolution of the shape 
of voids is a sensitive probe of the nature of dark energy. On the 
basis of a set of N-body simulations of structure formation in five 
different dark energy cosmologies, including dynamical dark en- 
ergy models, we confirm this sensitivity in the case of voids in the 
dark matter distribution. 

The first observation is that the SUGRA model has a less de- 
veloped void population and shows a lower level of clustering than 
the ACDM model, with the other quintessence models representing 
intermediate cases. The fact that the extended quintessence mod- 
els do not have a more manifest and pronounced signature of their 
rather distinct nature is one of our surprising findings. 

A key component of our assessment is in how far this finding 
is affected by the sparsity and spatial bias of the objects that probe 
the underlying dark matter distribution. In other words, whether 
this dark energy sensitivity may be exploited when studying the 
void population in the observed galaxy distribution. We find that 
the discrete, sparse and biased character of the halo and galaxy dis- 
tributions is seriously impeding the potential of probing the nature 
of dark energy from voids traced out by these objects. 



8.1 The a 8 degeneracy 

We have also looked deeper into the relation between the shapes of 
voids and the amplitude of the density fluctuations. Following the 
relations derived by Park & Lee ( 2007 ), we know that the ellipticity 
of voids is largely a reflection of a 8 , which quantifies the average 
amplitude of density fluctuations. In this study, we have shown that 
most of the effect is indeed degenerate with respect to cr 8 . 

In fact, we find that differences in the void shape (e)(z) be- 
tween the various dark energy models can almost be completely 
ascribed to differences in erg (z), i.e. the amplitude of density fluc- 
tuations. Of course, the differences in as and cosmic growth factor 
at any one cosmic epoch between the different cosmologies is the 
result of the differences in the equation of state parameter w(z). 
As is clearly borne out by figure [16] the amount of clustering fully 
determines the evolution of the void ellipticity. The fact that (e)(z) 
and as are so strongly correlated is interesting by itself. As they en- 
capsulate two different aspects of the large scale mass distribution, 
we would not immediately have suspected the existence of such a 
strong one-to-one connection. 

This leads us to the conclusion that if we wish to use void 
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shape as a probe of dark energy, we need to measure not only the 
void ellipticity (e)(z) but also, independently, cr&(z) at the same 
redshift. This will enable us to break the degeneracy between (e)(z) 
and 0-8(2). It ^ e a considerable challenge to improve the accu- 
racy of erg (z) measurements at all redshifts to anything compara- 
ble to that determined from the CMB. This also involves a bet- 
ter understanding of the non-linear effects on the evolution of ag 



( Juszkiewi cz et al.|2010) . Recent measurements of the growth fac- 
tor f(Ci) as a function of redshift are a great step in the right di- 



rection (e.g. 


Blake et aLpOU] 


Nusser, Branchini & Davis 2012, 


Tojeiro et al. 


2012| Blake et al. 


2012). 



8.2 Probing dark energy in the observational reality 

The discrete and sparse nature of the halo and galaxy distribution 
form a major source of confusion in retrieving the signal of void el- 
lipticity evolution. This adds to the confusion as a result of redshift 
distortions in the inferred galaxy distribution maps. 

We have tested the void ellipticities in a halo distribution 
whose average density is around 0.019/i 3 Mpc -3 , comparable to 
that of the SDSS DR7 galaxy sample ( Mo ntero-Dorta & Prada] 
2009). Volume limited samples will have even lower sampling den- 
sities. Sparse galaxy samples like these will render it very hard to 
extract any significant signal of the nature of dark energy. For the 
exploitation of void shape statistics in such observational circum- 
stances, significantly denser samples or samples over considerably 
large volumes will be needed. 

A second major complicating factor is the spatial bias of the 
halo and galaxy distribution with respect to the underlying dark 
matter field. Even when the observational data sets would consist 
of a densely sampled halo or galaxy distribution, they still may not 
form a fair reflection of the underlying mass distribution. Spatial 
bias is a major complication in the case of the halo population in- 
vestigated in this study. 

The situation would be even more complicated for the galaxy 
distribution. Baryonic processes involved in the formation of galax- 
ies in and around voids produce effects which are not yet fully 
understood. Recent work has shown that different semi-analytical 
prescriptions of galaxy formation (SAM) lead to significantly dif- 
ferent galaxy populations in and around voids (De Lucia etal. [ 2006 ; 
Bower et al. 2006). The same mass distribution may imply a biased 
void galaxy population in one SAM, while another implies an anti- 
biased population (Platen 2009 ; Platen, van de Weygaert & Jones 
2012). A possible solution to problems related to the spatial bias of 
haloes and/or galaxies would be to infer information from differ- 
ent unbiased sources. A promising possibility might be the use of 
(dark) matter maps determined by gravitational lensing. 

Another - practical - factor that influences our results are the 
instruments that we use for identifying and studying the void pop- 
ulations. The watershed void finder (WVF) succeeds beautifully in 
delineating the often erratic outline of voids. However, the DTFE 
density field reconstruction - from which the WVF voids are de- 
termined - may be quite noisy in the case of the sparsely sampled 
data encountered in the observational reality. We are in the process 
of investigating other density estimators and void finders. Prelimi- 
nary results show that the MAK reconstruction formalism used by 
|Lavaux & Wandelt ( 2010) may partially alleviate this practical con- 
sideration (Ruwen 2011). 

Circumventing the complications with spatially resolving the 
void population with the sparse galaxy or halo population is the 
suggestion by Lavaux & Wandelt ( 2011) to combine the signal of 
all sampled and observed voids via a scaled stacking of the voids. 



The resulting elongation or flattening of the stacked void may then 
be ascribed to the Alcock-Paczynski effect (Alcock & Paczynski 
|1979| ). The measured size differences in radial and transverse di- 
rection of the stacked void can then be directly related to the prod- 
uct of angular diameter distance and Hubble parameter. This deals 
to a large extent with the discreteness and sparseness of the data 
sample, and |Lavaux~ & Wandelt ( 201 1) argue and demonstrate con- 
vincingly that it leads to an highly accurate assessment of the dark 
energy equation of state (also see Shoji & Lee 2012). Nonetheless, 
it may still suffer from uncertainties on the biasing properties of the 
void galaxy population which, as we have seen in this study, may 
have a sizeable influence. 
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APPENDIX A: MODELS OF DARK ENERGY 

In this section, we elaborate on the specific models of dark energy 
that we have considered in this work (the same models were used 
in |De Boni et al.|(2011| >). Dark energy has its influence on cosmol- 
ogy through the Friedmann equation (equation^. Hence, we derive 
w*de(cl) for each model, which is the only remaining missing piece 
from equation [2] 

Our reference model is a universe containing cold DM and 
a cosmological constant: ACDM. We compare this model to four 
different models of time dependent dark energy. We use two 
quintessence models, in which the dark energy is described as a 
scalar field under the influence of a potential. The other two models 
are extended quintessence models, where the scalar field is coupled 
to gravity. 

In the following we set ao — 1. We assume a universe with flat 
geometry, i.e. without curvature. The equations used to determine 
w(a) are given. The resulting relations for the different models are 
shown in figure [2] 

Cosmological constant Dark energy in a ACDM cosmology is 
modelled by a cosmological constant, or equivalently a constant 
wa = —1 in equation [2] 

Quintessence Dark energy modelled by a scalar field in a po- 



Model 



ACDM 



RP 



SUGRA 



EQp 



EQn 



tential V(<j)) is called "quintessence" dark energy ( |Wetterich|1988 
Ratra & Peebles 1988}. This model has w = w(a) and the Fried 
mann equation is 

2 



u_ 



+Q o ,0exp -3 



Ja «' 



where 



Wcf, 



1^2 
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4> 2 - v{<j>) 
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Note that when the kinetic term vanishes, we regain the ACDM 
value of w = —1. The cosmological constant can, thus, be seen 
as a special case of the more general quintessence model of dark 
energy. We can solve for (j) using the Klein-Gordon equation: 

dV(</>) 



> + 3H(t> + 



0. 



(A3) 



The potential V(<j)) determines the model's dynamical proper- 
ties. We have used an inverse power law potential ( Ratra & Peebles 



1988 ) and a generalised inverse power law potential (Brax & Mar- 
tin 2000 ). The latter potential expands upon the former by includ- 
ing corrections from supergravity ( Freedman, van Nieuwenhuizen 
& Ferrara 1976 ). These models, which we will later refer to as RP 
and SUGRA respectively, have the following potentials: 

A 4+a 

VrpW = -t- (A4) 



w 



0.0 
-1.0 



0.0564 
-0.9 



0.452 
-0.9 



0.0117 0.0805 
-0.9 -0.9 



Table Al. Fits of dark energy model parameters w a to the w(a) relations 
in figure [5] determined using a x 2 fit. 



Extended Quintessence Secondly, we consider a scalar field ex- 
plicitly coupled to the rest of the universal components through 
gravity (Boisseau et al. 2000). Specifically, we consider here the 
so-called "extended quintessence" (EQ) models ( Pettor ino & Bac-| 
cigalupi 2008 ). The way we represent an interaction in field theory 
is by adding an interaction term to the action of the field. This term 
is a (Lorentz invariant) product of the quantities that represent the 
fields that we want to interact. In our case, these are the gravita- 
tional field represented by the Ricci scalar R and the EQ field (j). 
The action then becomes (Baccigalupi, Matarrese & Perrotta 2000 ) 



-F(<P)R--d^d^-V(<P) + L' 



where L' contains the terms of the Lagrangian without ( 
given by 

1 
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F{(j)) is 



(A7) 



where £ determines the strength of the interaction and 0o = (j)(to). 

In what follows, we only consider the linear (Newtonian) 
limit. We can then approximate the effect of EQ in an N-body simu- 
lation by replacing the gravitational constant G by a time dependent 
parameter G, given by: 



G 
G 



1 - SnG^ 2 ~ 4>l) • 



(A8) 



This is supported by version 3 of the GADGET code and, thus, 
we conveniently solve this otherwise complicated problem. Addi- 
tionally, in the linear regime, the equation of state parameter w(z) 
behaves like the normal quintessence ones. We use the RP potential 
for these models. 

The linear approximation is valid if wjbd ^> 1, where 



WJBD 



F{4>) 



4£ 2 



(A9) 



[dF{$)/d4>] 2 

Using this relation we can determine the allowed values of £. The 
lower limit for wjbd,o (and, thus, the upper limit for the interac- 
tion term £, because wjbd,o oc £ -2 ) can be determined by ob- 
servations. On cosmological scales, this limit - as obtained using 
WMAP1 and 2dF data - is set at \w JB d,o\ > 120 ( |Acquaviva 
et al. 2005 ). Because we want the interaction to be as strong as pos- 
sible within observational limits we indeed set £ using this value; 
wjbd,o — 120. 

The two models we consider are those with a positive and a 
negative value of £ (referred to as EQp and EQn hence forward). 
These two models differ slightly in the gravitational parameter. At 
z < 0, § > 1 for EQp and ^ < 1 for EQn. These corrections are 
within the few percent level. 



Vsugra(</>) = — exp (AtvGc/) 2 



(A5) 



where a ^ and A are free parameters. They are both tracker 
potentials. 



General parametrization We can fit equation [T] to these models 
using a simple % 2 fitting procedure with w a as a free parameter (wq 
is fixed to the values that were chosen for the simulations). We find 
the best fits as given in table [Al] (and again in table [TJ. 
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